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Abstract 

We consider the decay of a metastable domain wall. The transition proceeds 
through quantum tunneling, and we calculate in arbitrary number of dimensions the 
preexponential factor multiplying the leading semiclassical exponential expression for 
the rate of the process. We find that the effect of the motion in transverse directions 
reduces to a renormalization of the tension of the edge of the wall in the semiclassi- 
cal exponent. This behavior is similar to the one previously found for breaking of a 
metastable string. However this simple property is generally lost for spontaneous decay 
of higher-dimensional branes. 



1 Introduction 



Metastable domain wall solutions arise in models with spontaneously broken approximate 
symmetry. The existence of such solutions was shown from different points of view, for 
example, in [1, 2, 3]. The origin of a metastable wall can be illustrated in a model with 
potential shown in Fig.l, it corresponds to the interpolation between the same vacuum state 
at two spatial infinities, e.g. at z — — oo and z = +oo with the field winding around the 
'peg' in the potential. 

Such configuration is classically stable in the sense that the solution with certain topolog- 
ical number (winding number) can not evolve classically to a solution with different topolog- 
ical number. However, such a transition can proceed due to temperature fluctuation, when 
the path, corresponding to the solution can be lifted over the barrier or due to quantum 
tunneling (the motion in a classically forbidden region), hence the domain wall can undergo 
a decay at certain conditions. In this paper we restrict ourselves to the later case. 

The decay of a domain wall is analogous to decay of metastable vacuum[4, 5, 6] in 2 + 1 
dimensions. Indeed, if a hole with area A is created in a wall, the gain in the energy is 
e A. The barrier, on the other side, that inhibits the process is created by the energy /li, 
with £ being the perimeter of the hole and ji being a tension associated with the interface. 
Thus, the 'area' energy gain exceeds the barrier energy only starting from a critical size 
of the hole created, i.e. starting with a round hole of radius R = 2/i/e. Once a critical 
piece has nucleated due to tunneling, it expands, converting the domain wall. Therefore the 
probability of the transition is given by the rate of nucleation of the critical holes in the 
domain wall. The important difference, addressed in this paper, from the 2 + 1 dimensional 
false vacuum decay is that the domain wall can move in the transverse direction(s) and this 
motion affects the nucleation rate at the level of the preexponential factor. 

Recently we have considered a similar problem of calculating the rate of a transition 
between two states of a string with different tension [7]. It was shown that in a case of 
complete breaking (transition into nothing) of a string with tension e the essential effect 
of the motion in transverse dimensions reduces to a renormalization of mass parameter \i 
associated with the interface between two states, such that the dependence of the rate on 
string's tension has the same form as for two dimensions, where it coincides with well-known 
one for false vacuum decay or Schwinger process of charged pair production in external 
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electromagnetic field [8], expressed in terms of renormalized parameter hr 

As a result of the calculation to be presented in this paper we find that a similar behavior 
also applies to the decay of a metastable wall. Namely the essential effect of the motion 
of the wall in the transverse dimension(s) reduces to a renormalization of the boundary 
tension parameter /i in the expression for the false vacuum decay in 2 + 1 dimensions. It 
can be noted that such a reduction is not trivial and holds only for the transitions of strings 
and two-dimensional walls. We have explicitly verified that such behavior is lost in similar 
transition of higher dimensional branes, where the transverse motion produces effects that 
are not reduced to a renormalization of the tension of the interface. 

The decay rate of a QCD domain wall with tension e with further creation of the interface 
with tension \x was considered in Ref. [9] by adapting the expression for metastable vacuum 
decay in 2 + 1 dimensions: 



dr ( leu /i 3 



vzM-^T I (2) 



where V is the preexponential factor. This factor is found from a calculation of the path 
integral over small deviations from the semiclassical tunneling trajectory. In the limit of 
small £ the result of such calculation in a (2+1) dimensional theory, can be readily copied 
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from the corresponding expressions in the equivalent (2+1) dimensional problem of false 
vacuum decay [6] 

C 

V d =2+i = ^73, (3) 

with C being a constant independent of wall tension e, which can not be found from effective 
description of the theory. In other words, the constant C is determined by details of the 
underlying 'microscopic' model, but the dependence of the factor (3) on e is universal. 

In this paper we calculate the domain wall decay rate per unit area in arbitrary number 
of dimensions d, and find the result in a similar form, 

where again the constant C does not depend on e, and /ir is the tension of the interface that 
includes the renormalization effect of the wall fluctuations in the transverse dimensions. 

The material in this paper is organized as follows. In Sec. 2 we formulate the problem 
in terms of the effective Euclidean-space action, and in Sec. 3 we consider the separation 
of variables in the relevant path integrals with this action and also find formal expressions 
for these integrals in terms of products. A Pauli-Villars regularization for the products is 
introduced in Sec. 4, and the actual calculation is done in Sec. 5. In Sec. 6 we demonstrate 
that the dependence of the result on the regulator parameter is completely absorbed into 
the renormalization of the interface tension \i. Finally in Sec. 7 we present our results and 
conclusions. 



2 Euclidean action 

The low energy effective action for the problem at hand is given by Nambu-Goto action for 
two and tree dimensional objects 

S = fiA + eV, (5) 

with V being the world volume of the wall, while A is the world area of interface. Nontrivial 
classical solution (bounce), which defines the exponential behavior of the rate is empty sphere 
with radius 

R= 2 A (6) 

£ 

surrounded by the metastable phase (Fig.2). The action (5) is the low energy effective 
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Figure 2: Bounce configuration 

expression, and it does not take into account the thickness of the wall, or of the interface. 
Hence it is valid only while one can neglect the structure of the objects and consider them as 
having zero thickness. If the thickness of the wall is of order r$ then the natural mass scale 
associated with it is M = l/r . Therefore one can write the conditions of applicability of 
the action (5) 



k < M , £ > r , (7) 

with k and I being any momentum and length scales in the problem. For instance the radius 
of the bubble should be much larger then the thickness of the wall 

M R = ^ » 1. (8) 

The probability rate of the transition is given by the imaginary part of the ratio of the 
partition functions calculated around the bounce and the trivial configurations 

^ = ^Im^. (9) 
dA AT Zx K J 

The imaginary part of Z 12 arises from one negative mode at the bounce configuration. 
Furthermore, due to three translational zero modes, the numerator in Eq.(9) is proportional 
to the total world volume AT occupied by the wall, so that the finite quantity is the transition 
probability per unit time (the rate) and per unit area of the wall. 

In order to evaluate the relevant path integrals with the pre-exponential accuracy we use 
the generalized cylindrical coordinates, with r, 9 and ip being the spherical variables in the 
(t,x,y) space (of the bounce), and z being the transverse coordinate. We consider only one 
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transverse coordinate, since the effect of each of the extra dimensions factorizes, so that the 
corresponding generalization is straightforward. We further assume, for definiteness, that the 
space-time boundary in the (t, x, y) space is a sphere of large radius L, where the boundary 
condition for the wall is zir = L) = 0. The small deviations of the wall configuration from 
the bounce can be parametrized by the variation of the bulk: z(r,0,tp), the radial (/) and 
the transverse (£) shifts of the boundary: 

r(9, V ) = R + f(9, ip) z(R, 9, tp) = ((9, <p) , (10) 

In terms of these variables the action (5) can be written in the quadratic approximation in 
the deviations from the bounce as 

a 4lT r3 i 16 7T/! 3 fj, 

t 



- • + ^/dfi(AA/-2/ 2 + A(^) + 



£ - f r 2 drdQg ij d t zd 3 z , (11) 



where the tensor is the Euclidean metric tensor in spherical coordinates = r,9,ip), 
while g a b is induced metric on the sphere (a, b = 9,ip), dfl = sin 9 d9 dip. 

The action around a trivial configuration in the quadratic approximation takes the form 

4-7T f pL 

S 1 = —eL 3 + - r 2 dr dtt g^ d lZ d jZ , (12) 



3 2 Jo 

3 Separating variables in path integral 

One can immediately see that the variable corresponding to the longitudinal variations / of 
the boundary is decoupled from other variables. This implies that the path integral over / 
can be considered independently of the integration over other variables and that it enters as 
a factor in Z 12 . On the other hand it is this integral that provides the imaginary part to the 
partition function, and it is also proportional to the total space-time area AT. Moreover, 
this path integral is identical to the one entering the problem of false vacuum decay in (2+1) 
dimensions and we can directly apply the result of that calculation [6]: 

-L Im J Vf exp [f J dQ (g ab d a fd b f - 2 / 2 )] = , (13) 

where C is independent of e and it depends on parameters of underlying theory. The ex- 
pression for the transition rate thus can be written in the form 



dF C 16vr/i 3 \ Z 



exp ( -^f - ) , (14) 



dA £ 7 /3 3£ 2 Zl 
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with the path integral Zyi running only over the transverse variables ( and z 

Z12 = I V(V Zl Vz 2 exp (-S 12 ) 
and involving only the quadratic in these variables part of the action (11) 



(15) 



S 



12 



-e L + 



16 7T/Z 3 



+ 



f JdQg ab d a (d b C + £ - 



R 



r 2 drdQg ij dizdjZ . (16) 



3""" ' 3e 2 

In the same quadratic approximation the partition function Z 1 for the trivial configuration 
is given by 

Z 1 = Jvz exp (Si) (17) 

with Si given by Eq.(12) and the integral running over all the functions vanishing at the 
space-time boundary: z(L, 9, if) = 0. 

It was shown in [7] that the calculation of the ratio of this type of the partition functions 
can be reduced to a calculation of partition functions associated with boundary only. It 
is clear what is meant by the boundary partition function for bounce configuration. It 
is possible to define a similar object for the flat wall configuration in the following way. 
Although there is no relation between the flat wall configuration and the sphere with radius 
R = -jr, one can calculate the partition function Z\ by first fixing the transverse variable z 
at r = R: z(R, 9, <p) = ((9, ip) and separating the integration over the bulk variables, hence 
introducing by hands the boundary for trivial configuration. As a result on gets the ratio in 
the form 

-2l2(boundary) 



'12 



Z\ Z\ (boundary) 

with boundary partition functions given by 



£l2(boundary) = J T>( exp ~ J dtt ( A {2) ( + J dQ(d r Z, 



r=R 



and 



Z\ (boundary) = J £>C ex P J ^ C 3 



r=R 



(18) 



(19) 



(20) 



with the function z c (r,9,ip) satisfying the Laplace equation Az c = with the boundary 
conditions 

z c (R,6,<p)=C(9,<p), z c (r = L) = 0. (21) 

The operator is the angular part of the Laplace operator in 3d (the Laplace operator 
on a sphere). 
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One can find the complete set of these functions by expanding the boundary function 
((0,<p) in the series of spherical harmonics 

C(0,<p)=£4mlW0,¥O • (22) 
l,m 

In each of the partial waves these functions are then found as 

R l+1 

z c {r, 6, <p) = A lm — Y lm {d, <p) ■ (23) 

Substituting these functions to the equations (19) and (20) and performing integration over 
the amplitudes A[ m yields 

oo / 1 \(2i+l)/2 

^(boundary) = M g [ fini + 1)+eR{l + 1) ) ( 24 ) 

and 

oo / l \ (2i+l)/2 

^l(boundary) = ATg [ £i2(Z + 1) j • (25) 

4 Regular izat ion 

Clearly, each of the formal expressions (19) and (20) contains a divergent product, and their 
ratio is also ill defined, so that our calculation requires a regularization procedure that would 
cut off the contribution of harmonics with large /. In order to perform such regularization 
we use the standard Pauli-Villars method and introduce regulator fields Z a with the weight 
factors C a such that 

£C a = l, £0^ = 0, (26) 

a a 

for any n less then some finite number. The action corresponding to the quadratic part of 
the Nambu-Goto expression (5) for small Z a : 

S 12 = | J dn g ab d a ( a d b ( a + I f R r 2 dr dtl d t Z a d,Z a + M 2 a Z 2 a ) . (27) 

with M a being each regulator mass, which physically should be understood as satisfying the 
condition M a Mq and still being much larger than the relevant scale in the discussed 
problem, in particular M a R ^> 1. The regularized expression for the ratio of the boundary 
terms in Z 12 and Z 1 thus takes the form 

-i -l 



■ / 12(boundary) ^ ^ 



-2^1 (boundary) 



^-12(boundary) 
12(boundary) 



^1 (boundary) 
1 (boundary) . 



(28) 
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where we introduced the notation 7Z for the regularized ratio, and the regulator partition 



functions Z^l 



12(boundary) 



and Z 



(R) 

1 (boundary) 



are determined by the same expressions as in Eqs.(19) 
and (20) with function z c being replaced regulator functions counterparts Z ac which still 
satisfy the boundary conditions similar to (21): 



Z ac (R,9,ip) = ( a (9,(p) , 



(29) 



but are the solutions of the Helmholtz rather than the Laplace equation (A — M 2 )Z a = 0. 
The solutions of the Helmholtz equation fall off exponentially at the scale determined by 
M a , and for our purposes only the behavior near the sphere r = R is needed. For this reason 
we write the equation for the radial part of the l-th angular harmonic as 



Z' l ' + -Z' l - l -^-^Z l -M"Z l = 0, 



(30) 



where for a time being we have omitted the indices a and c for regulators. Introducing new 
function Z\. Zi(r) = J R/r Zi(r) we can rewrite previous equation in the form 



Z'l + -Z[- {l + y 2)2 Zi — M 2 Zi — 



(31) 



One can now write the radial coordinate as r = R + x, and treat the parameter (x/R) as 
small, since the scale for the variation of the solution is x ~ 1/^M 2 + (/ + l/2) 2 /R 2 . This 
approach yields an expansion of the regulator action associated with the boundary at r — R 
m powers of l/^(MR) 2 + (/ + 1/2) 2 . With the accuracy required in the present calculation, 
the (normalized to one at r — R) solution to Eq.(30) is found in the first order of expansion 
in (x/R) as 



Zi(R + x) 




1 - 



{MRf 



x 



2 (MR) 2 + (I + 1/2) 2 R 



exp -^(MR) 2 + (I + l/2) ; 



\x\ 



(32) 



Using the form of the solutions for the harmonics of the regulator field given by Eq.(32) and 
the expressions (19) and (20), one can write the regularized ratio of the boundary partition 
functions (28) as 



* = n 

1=0 



l(l+l) + 2 ] jM 2 R 2 +(l + ±) +1 
(/ + !)(/ + 2) 



(2/+1) C a /2 



X 
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n 



n 

1=0 



l + l 



(21+1) C a /2 



M 2 a R 2 +(l + \) +\_ 



1 + 



x 



M 2 a R 2 



M^ + (Z + I) 2 ) (l{l + l) + 2sjMlIV + (l + \) 2 + l 
where we took into account that R = 2 fi/e. 



(21+1) C a /2 



(33) 



5 Calculating the products 

Each of the products in Eq.(33) is finite at a finite M and can be calculated separately. 
Instead of calculating the product directly, we can use the relation 



n^z = exp(^>i^, 



(34) 



and calculate the sum. We start from the third product. The expression under the sign of 
product is of the form 1 + g(l), with g(l) close to for any /, since it behaves as M" 1 . Hence 
we can leave only the first term in the expansion of the logarithm ln(l + x) = x + 0(x 2 ). 
Thus, we need to find the sum 



=0 « 



M 2 a R 2 +(l + l) l(l + l) + 2^M 2 R 2 +(l + ^) +1 



(35) 

The sums associated with the three products are readily calculable with the help of the 
Euler-Maclaurin summation formula and the result for the regularized ratio has the form 



1Z = exp 



M 2 R 2 In MR + MR In MR + In MR 



(36) 



where M n lnM = ^ Q C Q M"mM Q , for any n. The expression for 1Z contains an essential 
dependence on the regulator mass parameter M. We will show, however, that all such 
dependence in the phase transition rate can be absorbed in renormalization of the parameter 
fj, in the leading semiclassical term. Although, it may appear that there is a problem with 
the term MR In M, which is not proportional to R 2 (the area of the interface) and thus is 
not the coefficient in front of tension //. But this behavior is only due to the particular 
dependence of the radius R on ji and e. 
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6 Renormalization of \i 



The parameter \i is defined in the action (5) as the coefficient in front of the area of the 
interface between the world space of the wall and empty space. Generally this parameter 
gets renormalized by the quantum corrections, and in order to find such renormalization at 
the level of first quantum corrections, one needs to perform the path integration using the 
quadratic part of the action around a configuration, in which the area of the boundary is an 
arbitrary parameter. For a practical calculation of this effect we consider a Euclidean space 
configuration, with the wall lying flat in (x, y) plane, and the interface being at y — 0. Thus, 
the area of the world surface of the boundary is XT, with T and X being the size of the 
world volume of the wall in the (t, x) plane. The Gaussian path integral over the transverse 
coordinates z(t,x,y) is then to be calculated. We use the notation ((t,x) fc> r the transverse 
shift of the boundary and expand it in the Fourier series 

C(M)=£a(£)e ife (37) 
k 

A similar expansion applies to the regulator boundary function £n(t,x). The part of the 
effective action associated with the boundary is determined by the functions z c (t,x,y) for 
the transverse shift of the string and the corresponding regulator functions Z c (t, x, y) that 
satisfy the equations 

Az c = and (A - Ml) Z c = , (38) 
and the boundary conditions 

z c (t, x, 0) = C(t, x) , Z c (t, x, 0) = ( R (t, x) (39) 

as well as 

z c (t, x, ±oo) = Z c (t, x, ±oo) = 0. (40) 
One can readily find these functions for each harmonic of the boundary values ( and Cr 

4(x,t,y) = e -^\y\j^ , (41) 

and 



Z*(x,t) =exp 



\y\ 



\A 2 + M* 



In order to separate the boundary effect in the path integral around the considered configu- 
ration from the bulk effects, we again divide it by the path integral around the configuration 
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where the whole world volume is occupied by the wall. Such division results, as previously, 
in the cancellation of the bulk contributions, and the remaining part of the effective action 
associated with the boundary is written in terms of the regularized path integral over the 
boundary function ( as 

fi R A = fiA- 

JV( exp{-| JdQ [fi g ab d a ( d b ( - e ( d y z c {x,t,y)\y^ +0 \} 



In 



In 



JT>Cr exp {-§ JdSl [ng ab d a ( R d b ( R - e(d y Z c (x,t,y)\y^ +0 \} 
JV( exp{-±eCdyZ c (x,t,y)\y^+o} 



fV( R exp{-±e( R dyZ c (x,t,y)\ y ^ +0 } ' 



(43) 



where /i R — /i + Sfj, is the renormalized mass parameter. The correction to ji can thus be 
written in the form 



5/i = 



d 2 k 



(2n) 2 
d 2 k 



1 

2 

I J (2tt) 2 

M 2 Me, e 2 

\ nM R + —\nMR + — — 

o7T 87T/i lb7TjJ, z 



In k 2 + -VM 2 + k 2 
. \ » 

In (M 2 + k 2 ) - In A; 2 



In k 2 + -Vk 2 



In MR. 



(44) 



It should be mentioned here, that result for the renormalized parameter /j, R would not change 
if we considered the wall not with a topology of a plane, but rather of a cylinder, namely, 
considering periodic boundary conditions in one direction. 



7 Results and conclusion 

Collecting all terms together we find the rate of the process. It is clear that the result for 
each d — 3 transverse dimensions factorizes, thus we have the expression for the rate in the 
form 

_ = _^-3 exp ^__fj, (45) 

where jj, is bare (non- renormalized) tension, and regularized ratio 1Z given by (36). Taking 
into account that each of the transverse dimensions contributes additively to 5/i and the 
interface is a sphere with area 

A = AttR 2 , (46) 
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and expressing the bare /i through the renormalized one: /i = fiR — <5/i, one readily finds that 
the dependence on the regulator mass M cancels in the transition rate, and one arrives at 
the formula given by Eq.(4). 

Thus we obtained the result similar to the decay of a string, when the effect of all extra 
transverse dimensions results only in the renormalization of parameter fi associated with the 
interface. It should be mentioned, however, that the result for a string was, actually, obtained 
for a transition between two states of a string with different tensions. Here we considered 
decay of a wall (transition into nothing). For the calculation used it is impossible to preserve 
finite terms, but only proportional to some power of the regulator mass parameter M. 

Having calculated the probability rate for a decay of one- and two- dimensional objects, 
e.g. string and domain wall, it is tempting to assume that the situation is somewhat similar 
for the decay of an object of arbitrary dimensionality. But it is not true already for the decay 
of tree- and four- dimensional objects, where the dependence of the result on regulator mass 
is substantial even after renormalization of a parameter associated with an interface. That 
dependence demands introduction of new terms into initial action, which corresponds to 
nonrenormalizibility of the effective 'low-energy' theory. 
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